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J3 '. ABSTRACT 

We present a method of cross-calibrating the polarization angle of a polarimeter using Bicep Galactic obser- 
vations. Bicep was a ground based experiment using an array of 49 pairs of polarization sensitive bolometers 

■*— > . observing from the geographic South Pole at 100 and 150 GHz. The Bicep polarimeter is calibrated to ±0.01 

in cross-polarization and less than ±0.7° in absolute polarization orientation. Bicep observed the temperature 
and polarization of the Galactic plane (R.A = 100° ~ 270° and Dec. = -67° ~ -48°). We show that the 
statistical error in the 100 GHz Bicep Galaxy map can constrain the polarization angle offset of Wmap W band 
to 0.6° ± 1.4°. The expected la errors on the polarization angle cross-calibration for Planck or EPIC are 1.3° 

■^j- ' and 0.3° at 100 and 150 GHz, respectively. We also discuss the expected improvement of the Bicep Galactic 

C . field observations with forthcoming Bicep2 and Keck observations. 
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1. INTRODUCTION 



The polarization of the cosmic microwave background radiation (CMB) provides a tool for probing the physics of 
the early Universe. The CMB polarization field is decomposed into even parity E-mode and odd parity B-modePJ 
Primordial density perturbations result in only E-mode polarization. The E-mode signal was discovered by 
Dasi and characterized by experiments including Boomerang, CBI, Maxipol, QUaD, Wmap, and BicepPEI 
Scientific interest in the CMB community moves toward detection of the B-mode signal, which originates from 
a primordial inflationary gravitational wave background and weak gravitational lensing. 3 Numerous kilo-pixel 
array experiments, including EBEX, Bicep2, Keck, PolarBear, Quiet, and Spider, are in operation or 
under construction to search for B-mode polarizationP3fl3E5EII 

While the sensitivity of an experiment increases by employing a large number of detectors, the requirements 
for controlling systematic effects becomes also stringent P^l Among systematic effects in the experiment, the 
polarization angle of the detectors is one of the most important quantities to be calibrated. Any miscalibration 
of the absolute polarization angle of a polarimeter mixes E-mode to B-mode signals, and therefore produces a 
false B-mode signal. Furthermore, such mixed E and B-mode signals are correlated and non-zero EB correlation 
indicates a false detection of the CPT violation or cosmic birefringencePSt251 
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The strategy for calibrating the absolute angle of a polarimeter may differ depending on the size of the tele- 
scope and the platform of the observatory, i.e. ground-based, balloon-borne and space-borne. The polarization 
angle of a small ground-based telescope like Bicep can be calibrated nearly end-to-end in the optical chain with- 
out replying on a calibration source on the sky but rather with a precisely oriented polarized source in front of 
the aperture in nominal observing conditions. On the other hand, large telescopes or any balloon- or space-borne 
telescopes are difficult to calibrate in nominal observing conditions without using a polarized sky signal. 

Commonly used polarized sources at millimeter wavelengths in the sky are the Crab nebula (Tau A) and 
Centaurus A (Cen A). The Crab nebula is a supernova remnant that emits highly polarized radiation. Aumont et 
al. presented the intensity and polarized signals of the Crab nebula observed by IRAM at 90 GHz.^U Cen A is an 
active galactic nucleus and Zemcov et al. reported the measurements of Cen A using the QUaD telescope.^ The 
reflection from the rim of the Moon is another source of the polarized calibration at millimeter wavelengths for a 
detector that has a large dynamic rangePS Wmap presented the measurements of the polarized celestial sources, 
including the Crab nebula, from 23 to 94 GHzpil The measurements of the Crab nebula show a consistent 
polarization angle with Aumont et al.. Planck satellite is planning to use this source to calibrate the polarization 
angle of LFI and HFI detectors P^ 

While these highly polarized compact sources are widely used for a polarization calibration, a high signal-to- 
noise diffuse Galactic polarized signal observed by Bicep is another polarized source for the angle calibration on 
the sky. Bicep was a millimeter- wave bolometric polarimeter that is designed to observe the CMB polarizationP3 
Bicep employs a refractive telescope with a small 24 cm aperture, simplifying the characterization of the end-to- 
end performance of the polarimctcr's entire optical chain. While the observation of Bicep is concentrated on the 
sky region that is minimally contaminated by the dust and synchrotron emissions, one-fifth of the observational 
time is dedicated to the Galactic plane observations. With systematic effects well controlled, a high signal-to- 
noise map of the diffuse polarized signal over the Galactic plane makes a standard calibration source on the sky 
for ongoing and forthcoming CMB polarization experiments. 

In Section [2] we discuss the statistical and systematic uncertainties in the Bicep polarized Galaxy map. In 
Section [3] we discuss the formalism to cross-calibrate the polarization maps produced by an unknown absolute 
polarization angle polarimeter using the Bicep polarization map. In Section [4] we apply this recipe to cross- 
calibrate between the Bicep 100 GHz and Wmap W band maps, as well as compute the expected constraint on 
the polarization offset angle for Planck and EPIC. 

2. BICEP POLARIZED GALAXY MAP 

Bicep was a ground based telescope observing from the geographic South Pole. The polarimeter consists of a 
two lens refractive telescope with a 24 cm aperture and 49 pairs of polarization sensitive bolometers (PSBs) at 
100 GHz and 150 GHz with a corresponding beam widths of 0.93° and 0.60°, respectively. A detailed description 
of the Bicep instrument is presented in Yoon et alp3 

Bicep observed two fields over the Galactic plane as shown in Figure [T] For each field, a telescope scans 
back-and-forth in azimuth at 2.8°/s over a 65° range at a constant elevation. The elevation is stepped by 0.25° 
after 50 right and left "half-scans" at a constant elevation. The telescope observed with four different orientations 
about its boresight: 0°, 135°, 180°, and 315°. Each observation of the single field has a fixed boresight angle and 
four observations cover all the boresight angles to increase the crosslinking coverage. 

2.1 Map making 

Figures [5] shows the Q and U maps observed by Bicep at 100 and 150 GHz. This section describes polarized map 
making, focusing on processes unique to the Galactic field analysis. The map making process that is common to 
the CMB analysis is described in Chiang et alP 

The low-level time stream cleaning is applied to the raw time stream in following steps, (i) deconvolution of 
a bolometer transfer function, (ii) low-pass filtering at 5 Hz, and (iii) downsampling to 10 Hz. The j th sample 
in a half-scan of a gain adjusted pair-differenced time stream is 

diti) = 1(*±M _ <M ), (i) 
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Figure 1. (a) The two Bicep Galactic regions are indicated over the FDS model 8 at 150 GHz! — ' Right: The integration 
time per nside=256 of Healpix pixel at 100 GHz (b) and 150 GHz (c)P 




Figure 2. The Q and U maps of the Galactic fields are shown in the unit of fj,K cm t,. Q and U are defined in the Galactic 
coordinate with IAU convention. 



where di A B are the individual low-level processed time streams of i th PSB pair, and g- lA B are the gain factors 
for each PSB calibrated at every elevation step. Common mode noise between the PSB pair, such as thermal 
fluctuations of the instrument and atmospheric fluctuations, is removed by differencing between the two PSB 
pair time streams, and we fit a third order polynomial, Fi(tj), to the pair-differenced time stream in order to 
remove residual 1// noise below ~ 0.1 Hz. 

When the telescope sweeps over the Galactic plane the variations in the time stream due to the Galactic signal 
and 1/ / noise are degenerate. To prevent removing the Galactic signal we apply a mask within the Galactic 
latitude of ±3°. The mask is applied in the time domain and we only fit the polynomial to the pair-differenced 



time stream where the mask is not applied. The httcd polynomial is subtracted from the time stream dij for all 
the samples in the half-scan. 

In some cases one end of the half-scan lies inside of the Galactic mask. We did not include such half-scan in 
the map making because the polynomial inside of the mask needs to be extrapolated from the edge of the mask 
and the extrapolated polynomial does not represent the 1// noise inside of the mask. The half-scans whose 
two ends lie in the mask are also excluded. Consequently, the recovered maps become "pac-man" shaped. The 
recovered region shrinks as the mask width increases. 

We follow the formalism described by Jones et al™ and used for the Bicep CMB analysis described in Chiang 
et alP The Q and U values of the pixel in the direction on the sky p are computed from d\, = di(tj) — Fi(tj) as 
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The weight lUy is an inverse of a variance of d\ ■ calculated from the samples outside of the mask in each half-scan. 
The angle ip is the PSB orientations projected on the sky and 7 = y^ is the polarization efficiency factor, where 
e is a cross-polarization response. Subscripts «a and «b refers to the A and B bolometers of the ith pair. We use 
nside=256 of Healpix pixelization to project Q and U on the skyP3 



2.2 Pixel noise in the map 

Statistical errors in the Bicep maps are estimated using jackknife maps. We split the data into three pairs of 
halves, (i) right and left half-scans, (ii) (0°,315°) and (135°, 180°) boresight angles, (hi) two sets of detector 
pairs located in alternating sectors of the 6-sector circular focal plane. We compute the Q and U maps of each 
data set as {mQ 11 mu 1 ) and (rnQ 2 ,mu 2 ). We compute the difference as 
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where N is the number of observations at each Healpix pixel. We compute the histogram of NEQ — Sjuq/ \/2f s 
and NEU = 5mu/V2f s (where f s = 10 Hz sampling rate) from the map pixels that meet the criteria of (1) 
N > 2000, (2) the Galactic latitude \6 g i a t\ < 3° and (3) not being at the edge of the observed regions. (Hereafter 
we call the region of the sky that meets these criteria a selected region.) We fit the histogram with a Gaussian 

Aexp 
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Figure [3] shows the histograms of the jackknife maps. The noise property is well described by the Gaussian 
distribution. The averaged NEQ and NEU from the three jackknives are 523 and 507 fiKy/s for 100 GHz, and 
428 and 431 fiKy^s for 150 GHz, respectively. Among the three jackknives the worst NEQ is 16 % larger than 
the best one that is from the right and left half-scan jackknife. 

2.3 Systematic error 

The polarization properties of a polarimeter are described by two quantities, PSB orientation ip and the cross- 
polarization response e. The PSB orientation tp is the angle at which the PSB is sensitive to the linear polarization. 
The cross-polarization response is the response of the PSB to the orthogonally polarized incident radiation. 
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Figure 3. The histograms of SrriQ and Smu for 100 and 150 GHz are shown for right/left (top), boresight angle (middle), 
and detector sets (bottom) jackknives. 





Figure 4. Comparison of the filtered Q (and U) values 
of the simulated map against to the input unfiltered 
Q (and U) values. The two linear lines are fits to the 
100 GHz (blue) and 150 GHz (green) data. The red 
line indicates the line that has a slope of 1 with zero 
offset. The offsets of 100 and 150 GHz Q signals are 
-5.6 and -18.1, respectively. 





Figure 5. The difference of the polarization angle at 
each pixel before and after applying the BICEP time 
domain filtering to the simulated maps is plotted. The 
points are selected from the pixels for the Galactic lat- 
itude of \6\ < 1 (red), \6\ < 2 and \0\ > 1 (black), 
|6>| < 3 and |6>| > 2 (blue). The four lines are [Aa| 
from Equation [9] with the cases for a?; n of 0.1°, 2°, 5° 
and 10°. 



Two calibration methods are used to measure the cross-polarization. One uses a modulated linearly polarized 
broadband noise source mounted 200 m away from the Bicep telescope. Bicep observed the source by raster 
scans with 18 different boresight rotation angles. The other method uses a rotating wire grid mounted at the 
cryostat window. The signal is generated by chopping between an ambient absorber and the sky. The cross- 
polarization response is measured to within ±0.01. 

The PSB orientation is measured with a rotating dielectric sheet mounted in front of the cryostat window 
in addition to the two methods described to measure the cross-polarization. The measurements were repeated 
through each observing year and the uncertainty of the individual PSB orientation is 0.1° rms. After the cryostat 
was opened between 2006 and 2007 observing years, the PSB orientation measurements showed an average of 1° 
rotation in the absolute polarization angle. Thus, the absolute PSB orientation uncertainty is assigned to be less 
than 0.7° rms for three years of the observation periods. The detailed discussion of the polarization calibration 
of the Bicep polarimeter is described in Takahashi et alP2 

2.4 Effects due to time domain filtering 

The subtraction of a polynomial fit from the pair-differenced time stream effectively acts as a high-pass filter in 
the time domain. While the purpose of the high-pass filtering is to remove 1// noise, this filtering also removes 
some modes of the signal. 



In order to quantify the amount of the time domain filtering of the Galactic signal, we prepare a simulated 
polarization map (Healpix pixelization of nside=256) that consists of the sum of CMB and FDS maps at 100 
and 150 GHz with a beam size of 0.93° and 0.60°, respectively.^ The CMB map is generated by synf ast using 
the cosmological parameters of the standard ACDM model presented in Komatsu et alpSEID The Q polarization 
of the FDS map is made based on the relationships of Q/T = co(T/T max ) Cl observed by BiCEPpS We have used 
(co,Ci) = (0.007,-0.47) and (0.017,-0.29) for 100 and 150 GHz, respectively. 

According to this polarization model, the Galactic Q depends on temperature signal. The FDS model 8 does 
not have the same level of the emission as it is observed by the BICEP. In order to simulate the realistic level 
of the Galactic emission we use the temperature T — fiTpDS, where j3 = 1.30 and 0.87 for 100 and 150 GHz, 
respectively. The U polarization of the FDS is set to be zero for all the pixels. The simulated maps are smoothed 
to the beam size of 0.93° and 0.60° for 100 and 150 GHz, respectively. 

We generate time ordered data using these simulated maps with the BICEP pointing and apply the same 
map making as we apply to the real data. Figure [4] shows the correlation between the input and filtered Q 
and U for the pixels inside of the selected sky region. The relationship of Q before and after applying the time 
domain filtering is well described by a simple linear relationship. The offset generally depends on the amount 
of the signal contained at the Galactic plane and the offset is higher when the Galactic signal is higher. This 
is because the signal level at the mask boundary is significant as compared to the 1// noise, and therefore the 
interpolated polynomial inside of the mask follows the trend of the Galactic signal instead of the trend from the 
1// noise. On the other hand, the U polarization does not show any clear trend. This is because U polarization 
do not contain any Galactic signal but only the polarization of the CMB. Therefore, there is no characteristic 
signal increase at the Galactic plane. 

Figure [5] shows the change of the polarization angle after time domain filtering. The change of the polarization 
angle |Aa| is modeled as 

Qfilt = Qin — Qo (7) 

UfiH = Ui n (8) 

Aq = -(arctan — arctan — — ) 

1 I p sm2a in 

= -(arctan — - a in ), (9) 
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where Qi n = I p cos2ai n , Ui n = I p sm2ai n and I p = \/Q\ n + Uf n . Qo is the offset to account for the filtering 
effect. The time domain filtering effect to the polarization angle ranges from 0.1 to 100 degrees. Therefore, when 
the offset angle between the Bicep map and the map from other experiment is cross-calibrated it is important 
to apply the same time domain filtering to the other map. 

3. ESTIMATION OF THE OFFSET ANGLE AND ITS ERROR 

We describe a method to detect the overall polarization angle offset between the two polarization maps. We 
have two sets of Q and U maps. Ones are the Bicep maps as the calibrated maps. The others are maps to be 
calibrated. In the case of comparing the maps from two different experiments, they do not necessarily have the 
same beam sizes, and therefore we need to deconvolve the original beam and smooth the two maps to the same 
beam size. The choice of the beam smoothing varies depending on the beam sizes of Bicep and other experiment 
to be calibrated. In this section we assume that the two sets of maps have a same beam size and are pixelized 
such that the noise among pixels are not correlated. We discuss the treatment of the different beam size between 
the separate experiments as a case-by-case basis in Section [4] 

We write the second and third components of the Stokes parameter of «th pixel in the two sets of maps as 

Bicep map : {Q iB ± 5Q tBl U lB ± SU iB ), (10) 

Uncalibrated map : (Qi ± SQi,Ui ± SUi), (11) 

where SQ and 6U indicate the statistical noise. We assume that the parent distribution of the pixel noise is a 
Gaussian described by the standard deviation of CTQ i(B) , <Ju i{B) , <Jqu U b) with a mean of zero. 



We relate Q and U of the same pixel on the sky between two experiments by two parameters, offset angle 5a 
and the ratio of the polarized amplitudes p as 
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We can solve Equation [T2l for pi and ai as, 
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and thus 
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where a B i = % arctan U B i/Q B i and a% — , arctan Ui/Qi. When the Q and U maps from two separate experiments 
are identical, a.i = and pi = 1. 

While the polarization calibration can be done in terms of Q and [/, we express Q and J7 of two maps in terms 
of da and p. This choice was made to mitigate the effect due to the spectral dependence of the instrumental 
bandpass location and shape mismatch between the two separate experiments We discuss the spectral dependence 
of the polarization angle in Section [S] 

When the Q and U maps contain only signals, we have a perfect knowledge of the offset angle 5a for every 
pixel. When the noise is present in the maps, the noise in the map has to be propagated to an error in the offset 
angle. The error of the offset angle in each pixel i is 
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where (TQ m , <7u B n a QU B n a Qn a Ui, OQUi ar e the pixel noise associated with Q B i and Ubi, and Qi and Ui, 
respectively. The derivative terms are 



ddai 
8Q B , 



1 



Um 



2Q|. 

1 



U%i 



Ui 



dSaj 



8Q t 2 Q 2 + U 



2 ' 



85 a.i 
dUs^ 
d5a 



1 



Q 



Bi 



2Q 2 Bl + 
1 Q t 



uh 



2Q +u 



(17) 
(18) 



The derivative terms are inversely proportional to the square of the polarized intensity. This indicates that the 
error of the offset angle is smaller when the polarized intensity is stronger. Figure [6] shows the angle uncertainty 
as a function of the pixel noise in Q and U maps and the polarized intensity, \j Q 2 + U 2 . 

While a polarization angle on varies from —90 to 90 degrees based on the signal and the pixel noise at the 
given point on the sky, we assume that the distribution of the differenced angle Sai is a Gaussian. We validate 
this assumption in Section |4] when we apply this formalism between Bicep and Wmap. 

The Galaxy is not a single point source, and therefore the estimation of the polarization offset angle improves 
by including all the available pixels in the map. In order to calculate the mean polarization offset angle, Sao, 
between the reference and the uncalibrated polarimeter maps and the corresponding uncertainty of the mean, 
we compute the likelihood of Sao a s 

Loce-^Vo), (19) 
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Figure 6. This plot shows the angle uncertainty of the 
polarization signal at a given pixel as a function of a 
pixel noise and a polarized intensity of the signal. For 
an example, a map that has a pixel noise of 10 \iK and 
polarized intensity of 100 [iK has a polarization angle 
error of 3° . 
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Figure 7. The spectra of BICEP (red) and WMAP 
W band (black). Both spectra are normalized to the 
maximum value of 1. The BICEP spectrum is the aver- 
age of a PSB pair. The WMAP spectrum is the average 
of the W band spectra. 



4. RESULTS 

We apply the method described in the previous section to Bicep and Wmap. 
constraint to Planck and EPIC, by using BICEP Galactic map. 



We also compute the expected 



4.1 Polarization angle offset between Bicep and Wmap 

Wmap has been observing the temperature and polarization over the full skypS The spectral bandwidth of 
Bicep 100 GHz and the Wmap W band overlap as shown in Figure [7J In this exercise, we assume that the 
absolute polarization angle of the Wmap polarimeter is unknown and we constrain the overall offset angle of the 
Wmap polarization maps using the Bicep Galactic map as a polarization calibration source. 

Before we apply the formalism described in Section |3l we need to correct the beam size difference between 
the two experiments. The FWHM of the Bicep beam size at 100 GHz is 0.93°. Each of the Q and U maps of the 
four Wmap W band differencing assembly is deconvolved with the corresponding Wmap P>i and convolved with 
FWHM of 0.93° Gaussian beam in nsidc=512 pixelization. We compute the weighted averaged map from the 
four differencing assembly maps in W band. The weights are the inverse of the pixel variance of each differencing 
assembly. We apply the Bicep time domain filtering to the averaged Wmap W band map. The filtered Q and 
U maps are downsampled to 0.92° pixel size (nside=64) maps in order to decorrelate the noise among pixels. 
The Bicep Q and U maps are also downsampled to the same pixelization. 

The pixel noise of the Wmap maps is computed by a — ,?° where <j§ for the Wmap W band differencing 
assemblies is (ow\ > <?w 3 > ^Wzi a Wi) — (5.940, 6.612, 6.983, 6.840) mK with nside=512 pixelization. We neglect the 
correlated noise between Q and U for both experiments. The pixel noise of the Bicep maps is computed based 
on the NEQ(U) derived from the right and left jackknife maps. 

Once the two sets of the maps and weights are computed in the same pixelization, we impose the criteria to 
select the pixels. We choose the pixels that meet the criteria of \0 g i a t\ < 3°, Nhits of Bicep > 2000 and pixels 
of which its neighbor do not have Nhits = 0. The second criterion assures that the most of the edge pixels of 
the map are not included. The second criterion does not exclude the pixel around 282° < 4> g i on < 322° and 
\8glat\ < 3° where the edge of the map is not tapered by Nhits- Therefore, we include the third criterion to 
exclude all the edge pixels in the maps. 

Figure [8] shows the map of offset angle Scti and the weight l/crf Q .. It is clear that the the offset angle is 
close to zero and the weight is higher at the Galactic plane. Figure [9] shows the weighted histogram of 8a%. The 
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Figure 8. The maps of Son (left) in unit of degrees and weight = 1/crf a . (right) in unit of degree 2 . The maps are 



downsampled to Healpix resolution of nside= 
not coincide with the ones in Figure [2] 



64. The edge pixels are removed, and therefore the shape of the map does 
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Figure 9. The weighted histogram of Sa between the 
polarization maps of Bicep and Wmap and the Gaus- 
sian fit are shown. The distribution is well described 
as a Gaussian. 
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Figure 10. The likelihood of Sa (solid black) with the 
pixel noise estimated from the right and left jackknife 
and Sa (dash black) with the pixel noise increased by 
16 % as a worst pixel noise estimation. The histogram 
is a mean of the Gaussian fit to da from the two sets of 
the simulated signal (CMB+FDS) at 100 GHz and the 
300 noise realizations. The red curve is the Gaussian 
fit to the histogram. 



mean and the standard deviation of the angle uncertainty of each pixel is —0.41° and 11.2°, respectively. The 
distribution of the histogram is well described as a Gaussian distribution. 

Figure [TU] shows the likelihood of the offset angle calculated based on Equation Q1J] using the Bicep and the 
filtered Wmap maps. The black line shows that the mean and the sigma are the 0.6° and 1.4° respectively. The 
dashed line with the same mean has 16% larger sigma as a worst case pixel noise. 

The histogram in Figure [TU] is the results of the signal and noise simulations. We prepare two sets of maps by 
adding the white noise of the Bicep 100 GHz and Wmap W band to the simulated signal only maps at 100 GHz 
described in Section 12.41 We repeat computing the mean of the Gaussian fit to the histogram of 6a from the 
two sets of the map for the 300 noise realization. The fit to this histogram in Figure \7\ is consistent with the 
likelihood obtained by the using Equation [TTJl 

4.2 Polarization offset angle between BICEP and future experiments 

Any ongoing and forthcoming CMB polarization experiments which observe the Bicep Galaxy region can cross- 
calibrate their polarization angle using the Bicep map. As examples, we compute the expected angle constraint 



Reference x Uncalibrated 


la error [°] 
100 GHz 150 GHz 


Bicep x No noise experiment 


1.24 


0.27 


Bicep x Wmap W-band 


1.45 


Bicep x Planck 
Bicep x EPIC-IM (4K option) 


1.26 
1.24 


0.27 
0.26 


(Bicep, Bicep2) x Planck 
(Bicep, Bicep2, Keck) x Planck 


1.26 
0.23 


0.08 
0.06 



Table 1. The la statistical error of the polarization angle offset for various combinations of the experiments. The first 
row, Bicep x No noise experiment, indicates the angle error only due to the Bicep statistical noise. Bicep2 only has 
150 GHz band, and therefore the error in 100 GHz does not show any improvement. 



for two cases, Bicep and Planck, and Bicep and EPICF^ 

Table 1 shows the list of la statistical error from the likelihood in Equation [i"9l for 100 GHz and 150 GHz for 
the two experiments. In this comparison, we assume that the bandpass shape of the two separate experiments 
is the same and the knowledge of the beam shape is perfect. 

The expected pixel noise of Planck and EPIC-IM are from Planck bluebook and Bock et al., respectivelyE3ESl 
It is clear that the estimate of the la error of the offset angle between Bicep and Wmap improves with Bicep 
and Planck or EPIC-IM. This is because the noise contribution from Planck and EPIC-IM is much smaller than 
the case from Wmap while the Bicep noise stays the same. On the other hand, there is negligible improvement 
from Planck to EPIC-IM because the source of the noise in these two cases is limited by the pixel noise of the 
Bicep map. 

While the observations of Bicep were completed, the ongoing Bicep2 and forthcoming Keck will improve 
the sensitivity to the angle calibration. If we assume that Bicep2 and Keck will spend the same observational 
time with the same detector sensitivity on the BICEP Galactic field, the expected reduction of the pixel noise is 

simply scaled by J N Btc 7 N , where Nsicep is the number of detectors of Bicep and Nq is of Bicep2 or Bicep2 
and Keck. We assume that N B ice P is 25 and 24 for 100 and 150 GHz, and N is 256 for Bicep2 150 GHz and 
144 x 4 and 256 x 2 for 100 and 150 GHz of Keck, respectively. The data combining with Bicep2 and Keck 
provide the statistical errors of the offset angle smaller than the systematic errors of the Bicep polarimeter itself 
for both 100 and 150 GHz bands. 

5. DISCUSSION 
5.1 Comparison between the diffuse Galactic source and the Crab nebula as a polarized 
calibration source 

We compare the Crab nebula and the BICEP Galactic region as a polarized source. The emission mechanism 
of the Crab nebula at the millimeter wavelength is dominated by the synchrotron emission. Macias-Perez et al. 
and Weiland et al. reported that the observed flux has a power law of ex ( 4n ^R, )~ ' 3 ~~ ' 35 while the degree of 
polarization stays constant around 7 % over the millimeter wavelengthPSEO On the other hand, the diffuse dust 
emission at the Galactic plane increases as a function of frequency.^ Therefore, the signal-to-noise increases as 
the bandpass location increases. 

The Crab nebula is a point-like source and the Galaxy is a diffuse source. In order to compare the two 
sources, we compute the integrated polarized flux of the Galactic field as shown in Table 2. We also show the 
integrated polarized flux reported in Weiland et al. and Aumont et alpSEl The spatial area of the Galactic 
field is much larger than that of the Crab nebula. Therefore, the integrated polarized flux of the Galactic signal 
within the BICEP field is larger than that of the Crab nebula. The total pixel noise from the Galactic field is 
added in quadrature. As a result, the polarized Galactic source at 150 GHz provides the same order of error as 
compared to the Crab nebula at W band. 



Band 


Q [Jy] U [Jy] Angle [°] 


WMAP, Crab nebula 
K band 
Ka band 
Q band 
Vband 
Wband 
all band combined 


-27.13 ±0.68 -1.40 ±0.08 -88.5 ±0.1 
-23.72±0.45 -1.88±0.12 -87.7±0.1 
-22.03 ±0.60 -2.06 ±0.14 -87.3 ±0.2 
-19.25 ±0.36 -1.52 ±0.24 -87.7 ±0.4 
-16.58 ±0.73 -0.75 ±0.42 -88.7 ±0.7 

0.07 


Aumont et al., Crab nebula 
90 GHz 


-88.8 ±0.2 


BICEP Galaxy 
100 GHz 
150 GHz 


114.7 ±6.0 29.1 ±6.2 ±1.5 
573.4 ±12.2 195.6 ±11.2 ±0.5 



Table 2. The integrated polarized flux of the Crab nebula and the Bicep Galactic region is shown. The polarization 
convention in this paper and Wmap are different, and thus the sign of U is changed from the original WMAP paper. 
The angle error of the Bicep Galactic measurements is the quadrature sum of the pixel noise. ^The polarization angle 
seen by 10' beam from Aumont et al. The original literature quoted the polarization angle in equatorial coordinates as 
a = 148.8°. 



5.2 Effect of spectral mismatch between two experiments 

When the two polarization maps from two separate experiments are cross-calibrated, the spectral bandpass of 
the two experiments is not necessary the same. We assess the effect of the bandpass mismatch to the offset angle 
estimation between the Bicep 100 GHz band and the Wmap W band. 

Gold et al. derived the synchrotron and dust emission templates by the Markov chain Monte Carlo fitting.^ 
We compute the simulated Bicep and WMAP maps by integrating the sum of the synchrotron and dust template 
maps over the Bicep 100 GHz bandpass and Wmap W band bandpass. We compute the offset angle 5a{ of each 
pixel between the two bandpass maps. The median offset angle of all the pixels within the selected sky region is 
0.005°. We define the signal-to-noise for each pixel as the ratio of the polarized intensity to the pixel noise. We 
compute the median and the maximum offset angle of which the pixels are the signal-to-noise > 3 are 0.01° and 
0.02°, respectively. 

The Bicep 100 GHz map expects a higher contribution of the dust emission as compared to the Wmap 
W band map because the Bicep 100 GHz bandwidth is slightly wider than Wmap W band bandwidth in higher 
frequency side as shown in Figure[7J Gold et al. shows that in the Bicep Galactic field the polarization direction 
of the synchrotron emission is —26° < a sync h < 0° and that of the dust emission is |«du.st | < 0.3°. Therefore, 
the overall offset angle between the Bicep and Wmap maps is expected to show the positive rotation due to 
the bandpass mismatch. The overall offset angle between the Bicep and Wmap maps, shown in Figure [7] , is 
0.6° ± 1.4°, and the positive mean value is consistent with the bandpass mismatch. 

This effect is prominent when the passband of the instrument is located where more than two emission spectra 
are mixed with nearly the same amplitude. This is because the two sources with different spectral shape can 
have different polarization angles. 

6. CONCLUSION 

We present the polarized diffuse Galactic emissions observed by Bicep at 100 and 150 GHz and the method to 
cross-calibrate the absolute angle between the Bicep map and any uncalibrated map. The absolute angle of the 
Bicep polarimeter is calibrated to ±0.7 ° and the la error of the polarization angle due to the pixel noise of the 
Bicep map is 1.24° and 0.27° for 100 and 150 GHz, respectively. 

We apply this method between the Bicep and Wmap W band maps and cross-calibrate the angle to 0.6±1.4°. 
The expected lc errors for the Planck 100 and 150 GHz bands are 1.26° and 0.27°, respectively. The ongoing 
and forthcoming Bicep2 and Keck are expected to reduce the statistical noise of the observations of the Bicep 
Galactic region significantly. 



The Bicep Galactic maps provide the polarized Galactic emission as a new angle calibration source for the 
ongoing and forthcoming CMB B-mode experiments that require the absolute angle calibration to a fraction of 
a degree. The method of using the Galactic signal as an angle calibration source can be applied to any two 
experiments if one of the polarimeters is well calibrated. Therefore, when the Planck full sky polarization maps 
are available, the future polarimeters should be able to use the Galactic signal as a calibration source not only 
with respect to Bicep but also to Planck. 
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